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ABSTRACT 

We address the study of the Ha vertical velocity field in a sample of four nearly 
face-on galaxies using long slit spectroscopy taken with the ISIS spectrograph 
attached to the WHT at the Roque de los Muchachos Observatory (Spain). The 
spatial structure of the velocity vertical component shows a radial corrugated 
pattern with spatial scales higher or within the order of one kiloparsec. The gas 
is mainly ionized by high-energy photons: only in some locations of NGC 278 and 
NGC 1058 is there some evidence of ionization by low-velocity shocks, which, in 
the case of NGC 278, could be due to minor mergers. The behaviour of the gas in 
the neighbourhood of the spiral arms fits, in the majority of the observed cases, 
with that predicted by the so-called hydraulic bore mechanism, where a thick 
magnetized disk encounters a spiral density perturbation. The results obtained 
show that it is difficult to explain the Ha large scale velocity field without the 
presence of a magnetized, thick galactic disk. Larger samples and spatial covering 
of the galaxy disks are needed to provide further insight into this problem. 

Key words: spectroscopy — corrugations— galaxies: spiral-kinematics. 


1 INTRODUCTION 

The Milky Way shows a galactic disk quite far from the 


planarity (Kerr et al. 19571. This high degree of struc¬ 


ture was first observed in the gas and can be classified 
into three categories associated with different regions and 
spatial scales of the disk: a) the gas in the innermost re¬ 
gion of the galaxy appears to be tilted with respect to 
the formal galactic plane (Gum et al. 19601 ; b) warp¬ 


ing of the gas in the outer regions of the disk is almost a 
ubiquitous feature of spiral galaxies in the Local universe. 


even for those that are seemingly isolated (i.e. Gum et al. 
1960 Verschuur||1973 Florido et al.||199~ l ; and c) since 
the 1950s we have known that HI in the Milky Way 
shows a coherent wavy structure with respect to the fun¬ 
damental plane of the Galaxy in the inner region of the 
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Sun called “corrugations” (i.e. Lockman 1977 Alfaro & 
Efremov|1996 1. 


Corrugations have also been observed in different 
stellar populations of the Milky Way (i.e. Quiroga|1977 


Alfaro, Cabrera-Cano, &: Delgado||1992 I and both in the 


azimuthal direction along different spiral arms (Spicker 
fc Feitzinger||19^l, and in the radial direction at differ¬ 


ent galactocentric radii ( Malhotra &; Rhoads|1995 1. This 
phenomenon is not just limited to our Galaxy but has 
also been observed in external galaxies (e.g. [Florido et al.| 


1991 Matthews & Uson 20081, and for different tracers 


of both gas and stellar populations. In addition, although 
most of the stellar corrugations were determined from the 
distribution of young stars, infrared observations seem to 
show that old stellar populations also participate in this 
structure ( jDjorgovski &: Sosin||198^ |Rhoads|1995| ). 

All these results suggest that this is a universal phe¬ 
nomenon common in late-type galaxies, observable at al- 
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most any wavelength (Alfaro fc Efremov|1996 Matthews 


& Uson 20081. However, despite the universality of this 


phenomenon it has barely been studied and is poorly 
understood. 

These spatial corrugations clearly must be associated 
with wavy vertical motions in the galactic plane. Evidence 
of these kinematic waves was first detected in the analy¬ 
sis of the rotation curves of spiral galaxies (e.g. |de Vau-] 
couleurs fc de Vaucouleurs|1963 Pi§mi§|1965 ), but it was 
not until Alfaro et al. ( 2001| analysed the velocity corru¬ 
gations in NGC 5427 in more detail that a likely physical 
mechanism was proposed for their origin. 

Henceforth, we use the term “velocity corrugation” 
to refer to the observed wavy field of the velocity ver¬ 
tical component. The two kinds of structures have been 
observed in spiral galaxies: a) spatial corrugations in our 
Galaxy and other edge-on external galaxies, and b) veloc¬ 
ity corrugations in face-on galaxies. 

Given the nature of the problem, galaxy-observer ge¬ 
ometry plays a fundamental role in the interpretation of 
the observables. Within our Galaxy we can easily ob¬ 
serve the corrugated structure of the spiral arms in the 
solar neighbourhood in the optical range, as well as the 
HI structure at larger distances in the radio wavelength 
range. The vertical velocity field of these morphologies 
shows an almost null component in the direction of the 


line of sight. However, Feitzinger & Spicker (19851 were 


able to detect true velocity gradients perpendicular to the 
Galactic plane that could not be totally explained by ge¬ 
ometric effects. 

By contrast, the face-on external galaxies represent 
the best natural laboratory for the study of the vertical 
velocity field, given that most of the observed radial ve¬ 
locity is representative of the galaxy disk velocity vertical 
component. However, they do not provide any informa¬ 
tion about spatial corrugations. 

In summary, in the last 60 years several pieces of evi¬ 
dence have been collected showing that most disks of spi¬ 
ral galaxies present a residual wavy structure with respect 
to the galactic main plane. This undulatory behaviour 
is observable at different wavelengths, being representa¬ 
tive of the gas in different physical conditions, as well as 
of different stellar populations. Different scenarios have 


been proposed to explain this phenomenon (see Alfaro fc 
Efremo\^|1996 Alfaro et akpOOl l, but there is a lack of 


systematic studies which would help to better understand 
these structures and their possible origin. 

In this paper, we address the study of the vertical 
velocity in a sample of four nearly face-on spiral galax¬ 
ies. As aforementioned, the geometry of the problem is a 
main criterion in the selection of the sample as well as 
of the long slit azimuthal position angles. In Section 
we show the galaxy sample and their selection criteria. 
Observations and data reduction details are presented in 


Section In section we describe the geometry of the 
problem and establish its fundamental relationships. The 
main results and relationships between the vertical com¬ 
ponent of the velocity field and the Ha emission lines are 
shown in Section]^ In Section|^we present an analysis of 
the ionization sources via Diagnostic Diagrams. Finally, 
a summary of this work can be found in Section]^ 


2 GALAXY SAMPLE 

The galaxies were selected from the LEDA databas^with 
the criteria: (i) low inclination angle (i< 25° (ii) nearby 
(v < 4000 km s“^); (hi) a well-defined spiral morpho¬ 
logical type (2<t<8); and (iv) a diameter larger than 2 
arcmin (D25 > 2). 

This information and other basic properties of the 
galaxy sample are summarized in Table With the aim 
of having a uniform and homogeneous criterion for the 
inclination and for the position angle, and since these val¬ 
ues can change significantly between different authors, we 


show the values given by Epinat et al. (20081. We also take 


the data for the calculation of the rotation curves from 
these authors. 

Figure [^displays images of the galaxies with the dif¬ 
ferent slit positions. Next, we summarize the main refer¬ 
ences and properties for these galaxies found in the liter¬ 
ature. 


• NGC 278 is isolated (Bushouse et al. 19881 and al 


most face-on, part of the local supercluster of galaxies 
centred around the Virgo cluster ( |Vallee|ri993p . It is rel¬ 
atively small in the optical emission, with size of about 
2' (Israel 2009) corresponding to a diameter of 7 kpc at a 
distance of ~ 12.1 Mpc (see Table [T|). However, i n HI its 
diameter is five times larger ( Knapen et al.|2d04 1. 

In optical images its morphology is that of a small, 
non-barred galaxy with a bright star-forming nucleus and 
spiral arms. Its nucleus is not only quite bright in Ha, but 
also in the infrared continuum (Dale et al.|2000 Sanders 


20031, and in the [CH] line (Malhotra et al. 20011, con¬ 


sistent with intense star-forming activity (Ho et al. 1995 
Eskridge et al. |2002[|. There is no indication of a low- 


luminosity AGN ( Tzanavaris fc Georgantopoulos|2007 |. 

NGC 278 can be found in many recent works, but 
mostly forming part of surveys (e.g. Israel [20TO I. A spe¬ 


cific study of this galaxy can be found in 


Knapen et al.| 


(20041, who concluded that its appearance most likely re¬ 


flects a recent minor merger. Asymmetries set up in the 


^ Lyon-Meudon Extragalactic Database, http://leda.univ- 
lyonl.fr/ 

^ For NGC 1058, LEDA gives i = 58.5°, but its correct in¬ 
clination angle is in the range 4° — 11° (Lewis 1987; van der 
Kruit & Shostak 1984; 2007AJ....134.1952P) 
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Figure 1. Images for the galaxies of the sample: NGC 278, NGC 1058, NGC 2500 and UGC 3574 ( HST images, WFPC2 
instrument and F450W filter). North is top and East to the left. The different slit positions of the long-slit spectra, taken with 
the red arm ISIS spectrograph attached to the 4.2m WHT, are also represented on the image . In the case of NGC 2500 an Ha 
image is also added, taken from the LVL survey, plotted as coloured iso-contour lines onto the HST image. The contours range 
from 2.5E-17 to 5.E-16 erg/s/cm^. 


disk by this minor merger may well be the cause of the 
intense star formation in the inner region, which can be 
interpreted as a rare example of a nuclear ring in a non- 
barred galaxy. Rather than being induced by the secular 
evolution of a bar, this nuclear ring would then be the 


direct result of an interaction event in the recent history 
of the galaxy. It also shows a disturbed velocity field in 
its Ha kinematics, which may be partly the result of spi¬ 
ral density wave streaming motions. The HI emission, al- 
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though faint, also shows disturbed morphology and kine¬ 
matics. 

• NGC 1058 is a nearly face-on spiral galaxy with an 


inclination angle in the range 4° to 11° (Petrie & Ru- 
pen|2007 1. The distance to this galaxy is uncertain, rang¬ 


ing from 8.4 Mpe (Battinelli et al. 20001 to 14.5 Mpe 
(Sandage fc Tammann||1974 (. We adopt a distance of 9.8 


Mpe, given by a more recent work (Moustakas & Kenni- 
cutt||200^|, for which 1" corresponds to 47.5 pc. 


It has a clear knotty structure not organized into spiral 
arms ( [Battinelli et al.||2000[ ), with a complete absence of 
any well-defined arms in Ha emission, and a fiocculent 


appearance of the disk’s innermost regions (Ferguson et 
jaTIP^ . These features are also observed in the HST 
image (Figure [^. 

Radial metallicity and reddening gradients, both de¬ 
creasing with the galactocentric distance, are found from 
studies of Ha emitting regions by [Ferguson et ^ (19981, 
a trend often observed in spirals ( Sanchez et al.|2014 (. 

The HI velocity dispersion (ranging 4-15 km s“^) also 
decreases with radius, but the decline of starlight with ra¬ 
dius is much steeper than that of the velocity dispersion 
( [Petrie fc Ruperi][2007[ ). These authors also find that the 
HI velocity dispersion is not correlated with star forma¬ 
tion or the spiral arms, so a global radial falloff must be 
explained in the context of significant local effects. 


Battinelli et al. {20001 also found a steeper radial slope 


for the star formation rate (SFR) in NGC 1058, with 
young star associations relatively small (peak at ~ 50 
pc), the star formation being clearly stronger in the cen¬ 
tral regions. However, from our long slit spectra. Figure 
Ha reveals stronger emission in the outer than in the 
innermost regions. 

• UGC 3574 is a nearby late-type spiral galaxy, with 


low inclination; at a distance of 21.8 Mpe (James et al. 


20041 it is the furthest galaxy of the sample. Two faint 


spiral arms can be distinguished, while the stronger emis¬ 
sion comes from the nuclear region and some very com¬ 
pact and clumpy sources at the end of the southern spiral 
arm or across the northern one. 

No specific studies for this galaxy are found in the lit¬ 
erature, but it is included in many surveys for different 
purposes, e.g. Boker et al. (20021 who study the nuclear 


star clusters in late-type spiral galaxies. In this work, 
UGC 3574 is characterized as having a “naked” nuclear 
cluster in the I-band WFPC2 image, forming a distinct 
and isolated entity within the disk and without any signs 
of a kinematic centre. 


• NGC 2500 belongs to a quartet of galaxies ( [Sandage] 
& Bedke|1994 1. Epinat et al. (20081 find diffuse Ha emis¬ 
sion in their observed Ha maps, in agreement with [James[ 
et al. (20041, where its short bar is almost aligned with its 


minor kinematical axis. There are no individual works for 


this galaxy, mainly found as part of surveys, e.g. [Epinat[ 


et al. (20081 or James et al. (20041. 

In the HST/WFPC2+F450W image. Figure]^ it is ob¬ 
served that the star formation is mainly concentrated in 
the nuclear region and its well-defined, northern spiral 
arm. The southern arm is fainter and, globally, the south¬ 
ern region of this galaxy appears more fiocculent. 

These features are enhanced in Ha emission, as can 
be observed in Figure The Ha image is taken from 
the public database of the Local Volume Legacy (LVL) 
survej]^ and is over-plotted to the HST image, in flux 
scaled contours. 


3 OBSERVATIONS AND DATA REDUCTION 

We obtained long-slit spectroscopy with the double arm 
ISIS spectrograph attached to the 4.2m William Herschel 
Telescope (WHT), at the Roque de los Muchachos Ob¬ 
servatory (La Palma) during December 2003. This in¬ 
strument consists of two intermediate dispersion spectro¬ 
graphs operating simultaneously, separating the blue and 
red spectral regions, which are imaged onto an EEV12 
and a MARGONI2 CCD chip, respectively. In this way we 
have two spectral ranges observed simultaneously, a blue 

o 

one centred around H/3 (4861 A), and a red one around 
Ha (6563 A). The gratings used, R1200R and R600B, 

o 

provide a dispersion of 0.23 and 0.45 A/pixel respectively. 
The slit width of 1 arcsec projects onto about 3.64 pixels 
Full-Width-Half-Maximum (FWHM) on the detector; the 
spatial sampling along the slit is 0.2 arcsec/pixel. The slit 
was placed at two or three different position angles for 
each galaxy (Table [^. An angle close to the minor axis of 
the galaxy was chosen, so the projection of the rotation 
velocity is nearly negligible and the observed velocities 
are mainly the projected vertical component on the line 
of sight. Moreover, other intermediate angles were taken 
to compare corrugated patterns with the latter, and/or 
to cross some brighter Ha features. 

The spectra were reduced and calibrated following 
the standard procedure. Bias subtraction, flat-fielding, 
wavelength and flux calibration were done with the IRAHj 
task ccdproc. For the wavelength calibration CuNe -I- 
CUAR lamps were used. The standard stars Feige 34, 
gl91b2b, gd248, gd50 and hz21 from the Oke (1990) cat¬ 
alogue were used for flux calibration. Sky subtraction was 
done using the IRAF NOAO package task background. 


^ http://irsa.ipac.caltech.edu/data/SPITZER/LVL/ 

IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agree¬ 
ment with the National Science Foundation 
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For each wavelength range, the different exposures at each 
slit orientation were combined, eliminating cosmic rays 
and bad pixels, thus obtaining the final spectra. 

The observed velocity along the slit is obtained from 
the Hq emission line fitting. Moreover, to calculate the di¬ 
agnostic diagrams to study the ionization (section]^, we 
fit the emission lines Ha, [Nil], [SII], H/3, and [OIII], with 
a Gaussian at each pixel along the slit where the emission 
was strong enough to be fitted. Line flux, FWHM, cen¬ 
tral wavelength, and their corresponding errors, are cal¬ 
culated with the STARLINI<0 package DIPSO. The con¬ 
tinuum level is simultaneously fitted locally with a first 
order polynomial. 


4 Vz CALCULATION 

In the case of face-on galaxies, the vertical component of 
the disk velocity field may be obtained straightforwardly. 
As our galaxies are nearly face-on, the observed velocity 
is the result of the sum of several projected velocity com¬ 
ponents, due to their inclination angle. This sum can be 
summarized as the vertical component to the disk, V±, 
plus its parallel component I^|. We assume then that the 
observed velocity can be expressed in terms of 


Vobs = Vsys + [I^i sin i + V± cos i\ = 

= Vays + [{Vrot COS 9 -f Vexp Siu 9) siu i + Vz COS i] (1) 

where "Vaya is the systemic velocity of the galaxy (radial 
velocity of the galactic centre with respect to LSR), Vrot 
and Waxp the rotational and expansion velocity respec¬ 
tively, at the plane of the disk, and Vz is the velocity 
perpendicular to the galaxy main plane; 9 is the angle in 
the plane of the galaxy (counterclockwise from the major 
axis), and i is the inclination angle of the galaxy disk. 

We consider the x-axis in the direction of the maxi¬ 
mum diameter passing through the galaxy centre and the 
north equatorial pole, contained in the plane of the sky. 
And the x-axis coincides with the major axis in the plane 
of the galactic disk. Thus, the polar coordinates in the 
plane of the sky, (r,(j>), and in the plane of the galaxy, 
(R,0), are related as follows: 

R = r Jcos2(0 - PA) + ~ (2) 

V cos i 

Rcos9 = rcos{(j> — PA) (3) 

Rsin9 = r sin(<() — PA)cos~^i (4) 

From equation neglecting the expansion velocity 


® http://star-www.rl.ac.uk/star/docs/sun50.htx/sun50.html 


(i.e., assuming Vexp ~ 0), and applying the relations be¬ 
tween the celestial coordinates (equations and 1^, the 
vertical velocity component can be calculated as: 

Vz = (Voba — Vaya “ ^ COs((/) — PA) siut^COS”^ i = 

( .. \ 

Vo6s "^sys Vrot “ 


V 


' L^ t^Gi4>-P ~^ 

V cos i / 


= (Vobs — Vsys) cos — Vrot cos 0 tani (5) 

where (from equations and 


cos 9 = — cos{4> — PA) = 

R 


cos((/) — PA) 

aj cos^ {<t> — PA) + sin^ [4’ ~ PA) / cos^ i 
= (l-f tan^(<() — PA)/cos^ i) 


( 6 ) 


The vertical velocities are therefore computed as a 
function of V ay a, the systemic velocity of the Galaxy taken 
from the literature; Vobs and (f>, the observed velocity, ob¬ 
tained by fitting the Ha emission lines from our long-slit 
spectra, and their corresponding slit position angles; and 
Vrot, PA and i, the rotational velocity of the galaxy, the 
position angle of the line of nodes and the inclination an¬ 
gle of the galactic disk. 

As a first approximation, we consider a global ve¬ 
locity model where the main movement in the galactic 
plane is the rotational velocity, neglecting any other kind 
of kinematics, and it is assumed that they depend on the 
galactocentric distance or radius, Vrot{R)- We try to ex¬ 
plain those movements that do not fit to the model as 
vertical motions with respect to the galactic plane. 

The observed velocities Vobs are obtained from the 
Ha emission line fitting of the long-slit spectra, described 
in Section 13 above. 

Data for the calculation of the rotation curves 
Vrot{R), inclination angles and PAs for each galaxy are 
taken from [Epinat et al.| ( |2008[ ), see Table Rotation 
curves data are given as a set of points with their errors 
(Figure , so the rotation curves are obtained by fitting 
these sets of points using the parametric model of |Gio-| 
|vanelli fc Haynes (20021, 


Vpe(r) = Vo(l - e ’'^’'*’')(1-I-ar/rpe) (7) 

where Vo regulates the overall amplitude of the rotation 
curve, Tpe yields a scale length for the inner steep rise, and 
a sets the slope of the slowly changing outer part. After a 
least square fit of this model to the data, the parameters 
obtained for each galaxy are summarized in Table 
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Table 1. Galaxy Parameters'^ 


Galaxy 

RA (J2000) 
h m s 

Dec. (J2000) 

o / n 

Type 

Redshift 

Dist.* 

(Mpc) 

Inclin.° 

(deg) 

PA= 

(deg) 

Dimensions 

(arcmin) 

Mfl 

NGC 278 

00 52 04.3 

-1-47 33 02 

SAB(rs)b 

0.002090 

12.1 

21±14 

52± 3 

2.1x2.0 

-19.6 

NGC 1058 

02 43 30.0 

-1-37 20 29 

SA(rs)c 

0.001728 

9.8 

6 ± 15 

125± 6 

3.0x2.8 

-18.7 

UGC 3574 

06 53 10.4 

-1-57 10 40 

SA(s)cd 

0.004807 

21.8 

19 ± 10 

99± 3 

4.2x3.6 

-18.0 

NGC 2500 

08 01 53.2 

-1-50 44 14 

SAB(rs)d 

0.001715 

11.0 

41± 10 

85± 2 

2.9x2.6 

-18.2 


Sourced from NASA Extragalactic Database 
^ Moustakas &; Kennicutt 2006; except for UGC 3574, James et al. 2004 

^ Epinat et al. 2008. Inclination and position angle deduced from their analysis of the velocity field (from Ho data). 

The galaxies were selected from the LEDA database^ with a criterium of a low inclination angle ( i< 25°). NGC 2500 
appears in this catalog with an inclination angle of 18°, but we use data from Epinat et al. (2008) to determine the galactic 
rotation curve, where a kinematic inclination angle of 41° is estimated. For the sake of coherence in our data, we take this 
value into our analysis. 

Table 2. Journal of observations 


Galaxy 

Date 

P.A. (deg) 

O 

Wavelength range (A) 

Exposure time (s) 

Slit (") 

Airmass 

NGC 278 

15.12.2003 

40.5, 93.15, 143.5 

6034-7088 

3x1200 

1.03 

1.06 - 1.33 

NGC 2500 

15.12.2003 

24, 79.5, 40.7 

6034-7088 

2,3x1200 

1.03 

1.08 - 1.14 

NGC 1058 

16.12.2003 

32.5, 125 

6034-7088 

3x1800 

1.03 

1.01 - 1.17 

UGC 3574 

16.12.2003 

0, 27, 122 

6034-7088 

3x1200; 2,3x1800 

1.03 

1.14 - 1.28 


Table 3. Parameters from the rotation curve fit for each galaxy. 


Galaxy 

Fo 

(km s“^) 

'h'pe 

(arcsec) 

a 

R 

(correlation) 

NGC 278 

219.73±59.1 

1.08±0.36 

-0.40±0.17 

0.907 

NGC 1058 

74.5±12.7 

0.16±0.10 

0.05±0.03 

0.668 

UGC 3574 

155.5±18.5 

1.29±0.31 

0.024±0.017 

0.916 

NGC 2500 

55.9±3 

0.40±0.53 

0.04±0.007 

0.872 


** The rotation curve has been fitted by a least squares fit, following the parametric model of|Giovanelli & Haynes| (|2002||: 
Vpe(r-) = Vb(l - + ar/rpe). 


Figure shows the rotation curve fitted to the data 
of each galaxy. The approaching side of the curve is rep¬ 
resented by grey circles and the receding side by black 
squares. Both sides of the curve were fitted separately, 
represented with dashed and dotted lines, respectively. 
The kinematic lopsidedness commonly observed in the dy¬ 
namics of galaxies can be clearly seen. This lopsidedness 
occurs in two different ways, as already described by |van| 
[Eymeren et al.| ( |2011[ ) , who found five different types for 
the kinematic behaviour of disk galaxies based on the ro¬ 
tation curves of the receding and approaching sides sepa¬ 
rately. NGC 278 shows their third type, in which receding 
and approaching sides agree well at small radii but differ 
at large radii, whereas the other three galaxies show type 


5, in which the curves of receding and approaching sides 
change sides at a certain radius. 

Despite this lopsidedness, the rotation curves were 
fitted to both receding and approaching points as a whole 
data set, given by the bold line in Figure This global 
fit remains at intermediate values between the rotation 
curves of the receding and approaching sides, as a combi¬ 
nation of both sides. 

Finally, as Vz is expressed as a function of 
and Vrot^SVrot (and PA±5PA are consider part of 
our data), we can calculate its uncertainty by propagating 
the error as: 


SVz = \cos ^i\SVobs + \—cos9taxi i\5Vrot (8) 
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Rotation Curve NGC 278 



Rotation Curve UGC 3574 



r(kpc) 


Rotation Curve NGC 1058 



Rotation Curve NGC 2500 



r(kpc) 


Figure 2. Rotation curves for the galaxy sample. Data are taken from |Epinat et ar] | |2008l l, where the approaching side of the curve 
is represented by grey circles, and the receding side by black squares. Both sides are fitted separately, represented with dotted and 
dashed lines respectively, and as a whole data set plotted with the bold line. This last fit is assumed as the rotation curve for the 
galaxy. Lopsidedness between the rotation curves of the receding and approaching sides is observed. 


Figurej^shows the procedure for obtaining the verti¬ 
cal motions for the galaxy NGC 278, at each of the three 
position angles. The top panels show the observed radial 
velocities Vobs with their error bars, together with the 
projection of the rotation curve Vrot (the bold line) fitted 
as described above. The systemic velocity Vaya is marked 
by a black-filled circle. All these plots have the galacto- 
centric distance in the a:-axis, where zero distance marks 
the galactic centre position. The Ha emission line flux is 
also plotted, with the aim of looking for some relation or 
correlation between the velocity and flux peaks. 

The middle panels show that Vz has some residual 
trend, as residual of the rotational velocity removal prob¬ 
ably due to the described lopsidedness. To remove these 
rotational residues, the vertical velocity component has 


been detrended, hereafter denoted as AVz, by fitting a 
linear component to the gentle rise for each rotation curve 
side, the approaching and the receding sides. In this sim¬ 
ple way, the residual trend from the difference between 
the approaching and the receding sides of the rotation 
curve is removed; and thus only the local fluctuations, or 
oscillations, of the vertical velocity component across the 
disk, throughout the spiral arms and for different slit po¬ 
sitions, remain. The bottom panels show the detrended 
velocities AVz. 

The same procedure is applied for the rest of the 
galaxies, although, for the sake of clarity, only the ob¬ 
served velocities (on the top panels) and the final de¬ 
trended vertical velocities (on the bottom) are shown in 
Figures]^ andIn the next section we describe a more 
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8 M.C. Sanchez-Gil et al. 




Figure 3. Top: The observed velocities for NGC 278, derived from the Ho; emission lines (see Section]^. The bold line corresponds 
to the projected rotation curve, and the black dot marks the systemic velocity. The Hu intensity is drawn in the right-hand ordinate, 
with a dashed line. The x-axis represents the distance to the reference pixel, in kpc units; Centre: The perpendicular velocity Vz, 
at each pixel position across the slit, calculated as Equation]^ Bottom: The detrended perpendicular velocity, by fitting a linear 
component to the global trend. A dotted line marking the zero value stands out as the global trend has been practically removed, 
with only the local oscillations of Vz remaining. The black line is a smoothing line of AV^, for an easier comparison with the Ha 
emission lines. © 2015 RAS, MNRAS 000, [T|jl9| 
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0 ) 


Figure 4. Top: The observed velocities for NGC 2500, derived from the Ha emission lines (see Section]^. The bold line corresponds 
to the projected rotation curve, and the black dot marks the systemic velocity. The Ha intensity is drawn in the right-hand ordinate, 
with a dashed line. The x-axis represents the distance to the reference pixel, in kpc units; Bottom: The detrended perpendicular 
velocity AV^, at each pixel position across the slit, calculated as Equation]^ The black line is a smoothing line of 
easier comparison with the Ha emission lines. 
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Figure 5. NGC 1058, PA = 32.5°. Top: The observed velocities for NGC 1058 and UGC 3574, derived from the Hq emission 
lines (see Section]^. The bold line corresponds with the projected rotation curve, and the black dot marks the systemic velocity. 
This PA slit position is very close to the minor axis and the rotation curve projection appears to be rather constant, equal to the 
Vsys, since the rotational velocity is negligible. The Ho intensity is drawn in the right-hand ordinate, with a dashed line. The 
x-axis represents the distance to the reference pixel, in kpc units; Bottom: The detrended perpendicular velocity by fitting a linear 

component to the global trend. The black line is a smoothing line of AVz, for an easier comparison with the Ha emission lineS; _ 
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detailed analysis or description of the morphologies or 
patterns found from these plots for each galaxy. 


5 RESULTS 

In this section we describe the deviations from the main 
rotation observed, AVz, which we interpret as perpendic¬ 
ular motions of the ionized gas from the galactic disk, or 
corrugations. The aim is to show the existence of such 
corrugations in galaxies, as well as giving a description of 
its morphologies or patterns. We found the maxima and 
minima values of AVz are usually related with Ha emis¬ 
sion line peaks, and with the morphology of the galaxy. 
Its interpretation is still an open matter of discussion. 

For NGC 278 and NGC 1058 a wavelike structure of 
AVz, along the whole slit, is clearer than for the other 
two galaxies. In NGC 2500 and UGC 3574 there are some 
regions with no Ha emission or too faint to be fitted, and 
consequently there is no AVz associated with these re¬ 
gions. This last feature is more remarkable in UGC 3574, 
the most distant galaxy of the sample, losing therefore a 
global view of the vertical velocity deviations with respect 
to the Ha emission (Figure]^. 

In a previous study for the corrugated velocities in 
NGC 5427, Alfaro et al.| ( |200l| found a systematic spatial 
shift between the velocity peaks, AVz, and the emission 
line intensity peaks of Ha; where the approaching (nega¬ 
tive values) peaks of AVz occur in the convex border of 
the arms, and the receding (positive values) maxima are 
located behind the Ha emission maxima, in the concave 
side. This kinematical behaviour is interpreted as the re¬ 
sponse of the gas flow into a spiral density wave in a thick 
and magnetized gaseous disk, described by the |Martos &:| 
Cox (19981 and Martos et al. (19991 galactic bore models. 


For each slit in all galaxies, we find some of these 
features of the variation of the vertical motions when the 
slit passes through spiral arms, but with a varied phe¬ 
nomenology. 

In NGC 278, the behaviour described by ([Alfaro et al.| 


20011 is clearly observed for the slit PA = 40.5°, between 
galactocentric distances —1 and —0.5 kpc, with AVz ~ 30 
km s“^ (Figure [^. This feature is in the concave side 
of the spiral arm, since it is located within the corota¬ 
tion radius (~ 1.5 kpc, Buta & Zhang 20091. The gas 


reaches the density wave behind the spiral arm border -or 
concave side-, so the approaching (negative values) peak 
in AVz occurs behind the Ha emission maxima, or the 
corresponding spiral arm, when the gas is flowing up to 
higher galactic latitudes above the arm. Then it falls down 
onto the galactic plane, which corresponds to the receding 
maxima (positive values) in AVz that appear ahead of 
the arm. 

This phenomenon is also observable, although for 


fainter Ha emission maxima, and close to the corotation 
radius, between galactocentric distances —1.2 and —1.5 
kpc, with high vertical AVz ~ 40 km s“^, and around 1 
kpc with AVz ~ 15 km s“^. 

Another interpretation for the approaching and re¬ 
ceding peaks in AVz could be the deceleration and accel¬ 
eration with respect to the rotation velocity. If this were 
the case, the observed velocities would be represented by 
Vrot- If the vertical motions were neglected, the associated 
variation in the observed velocities of AVobs ~ 20 km s“^ 
cannot be explained by the difference between the rota¬ 
tion velocities. 

An even more outstanding example is found for the 
slit PA = 143.5°, the minor axis. The most important 
peak in this case shows a variation in the observed veloc¬ 
ities between galactocentric distances 0.5 and 0.8 kpc. A 
AVobs ~ 22 km s“^ cannot be explained by the difference 
between the rotation velocities A V^ot cos 0 sin f ~ 0.13 
km s“^. On the other hand, if the AVz were not consid¬ 
ered and with a projection factor of cos 0 sin i ~ —0.009, 
for a AVobs of 22 km s“^ a difference for the rotation 
velocity of the order of 2400 km s“^ would be required. 
This would imply a velocity radial gradient ~ 8000 
km s“^/kpc, which makes no sense in a galactic disk and 
radial distances of 300 pc. 

With our data and resolution we cannot discern if 
this is the case of a bubble in expansion or a merger, as 
the process that would explain such behaviour in AVz. 
Moreover, this feature at PA = 143.5° is found as the 
largest vertical displacement, whereas it does not seem to 
be related with a significant Ha emission peak. However, 
we cannot claim that it is not related with a spiral arm, 
since the weakness of Ha emission does not determine the 
presence of a spiral arm. 

At slit PA= 93.15°, between ~ —1 and 1.2 kpc (close 
to corotation), we can identify a double AVz peak associ¬ 
ated with Ha emission peaks. So its interpretation in the 


framework of Alfaro et al. (20011 is not clear. Although 


this is a quite plausible model, with the current data, we 
cannot discard other mechanisms as the origin or con¬ 
tributing to the enhancement of the corrugations in the 
velocity field. 

Similar patterns can be observed in the rest of the 
galaxies. In some locations, it is possible to discern similar 
behaviour to that expected from a galactic bore. For many 
others, the mechanism behind it is not clear. However, the 
changes in AVz are still clearly seen, as it crosses a spiral 
arm (as identified by the Ha peaks), from the receding to 
the approaching side, or vice versa. So these changes in 
AVz are definitely related to the spiral arms. 

Some examples of such patterns are found, for exam¬ 
ple in NGC 2500, between —2.5 and —2 kpc at PA= 24°, 
AVz goes from —10 to 20 km s“^. Again, just after the 
spiral arm, the gas flows up to nearly 30 km s“^. At 
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PA= 40.7°, in the central part of the slit, AVz are more 
pronounced but they are also noisier due to the faint Ha. 
The stronger Ha emission peaks are at the ends of the 
slit, where we can observe the pattern described above, 
around —2 and ~ 1.8 kpc, respectively. At PA= 79.5°, 
close to the minor axis, this is observed between 3 and 
3.5 kpc, where the gas flows up from 30 to —20 km s“^. 
Between the two arms there is a similar change in AVz 
from —10 to 30 km s“^. As in the previous PA, there is 
a weaker and noisier zone between —2 and 1 kpc, but the 
corrugations are still observable. 

In NGC 1058, for PA= 32.5° and around the corota¬ 
tion radius at —1.5 kpc, there is another example of AVz 
across the spiral arm. In the opposite side of the slit, be¬ 
tween 1.5 and 2 kpc, just before the spiral arm there is an 
important change of AVz = —10 to 15 km s“^. Another 
example of the pattern we are describing as similar to a 
galactic bore is found at PA= 125°, between 1 and 1.5 
kpc, just before the arm (around 1 kpc) where there is 
a large vertical deviation. At this PA the Ha emission is 
bright enough that peaks have been measured continu¬ 
ously, and so corrugations are more noticeable. 

UGC 3574 is the farthest galaxy of the sample, there¬ 
fore the Ha emission is fainter and more difficult to mea¬ 
sure along the slit (Figure]^. Nevertheless, the vertical 
velocity deviations are appreciable and these seem to be 
linked to the spiral arms. At PA= 27°, between —9 to 
—8 kpc and 1 to 2 kpc we can see how the AVz changes 
across the spiral arms. 

In summary, we found in some cases that it is easy 
to see how AVz changes across the arm, quite similar to 


took the possible corrugation of velocities into considera¬ 


the pattern described by Alfaro et al. (20011, and it could 


be explained by galactic bore models. But in many other 
cases, AVz changes are shifted, or they become wider, 
or lengthened beyond the arm. That is, there are other 
features where the expected pattern for this model is not 
so clear or is not observed. In some cases, an important 
vertical displacement AVz does not seem to be related 
with an Ha emission peak, as for example the feature 
described above at PA = 143.5° in NGC 278 or between 
— 1 and 0 kpc at PA = 40.7° in NGC 2500. 


tion was published by Alfaro et al. (20011, for the galaxy 
NGC 5427, in which the hydraulic bore mechanism ( |Mar-| 
tos fc Cox|1998 Gomez fc Cox|2004a|b I was proposed for 
the formation of the velocity field observed in the vicinity 
of the spiral arms. The choice of this model, for compar¬ 
ison with the observations of the velocity field of NGC 
5427, was determined by it being the only model that 
predicted the expected velocity structure in the neigh¬ 
bourhood of spiral arms, and therefore the only one that 
made it possible to perform the comparison. Moreover, 
both the quantitative and qualitative aspects of the model 
fitted with the velocity field observed. The situation has 
not changed much over the last decade, and few models of 
the Galaxy, if any, incorporate corrugations into their de¬ 
scription. Yet, what is the principal ingredient that makes 
the models proposed by the Wisconsin group (and suc¬ 
cessors) different to the rest? The answer is simple: the 
magnetic field. The majority of models prior to the one 


proposed by Martos & Cox (19981 are based on the in¬ 
teraction between different galactic subsystems (like discs 
and halos), and the impact of the disc with high-velocity 
clouds (whether HI or dark matter), and do not include 
the magnetic field in their physical fundamentals - with 
the exception of Santillan et al. ( 1999| ) on the impact with 
HVCs, and those on the Parker instability for the gen¬ 
eration of super-clouds and corrugations in spiral arms 
(i.e. Franco et al. [20(121. Should we include the magnetic 


field in the physical foundations of any model that tries 
to explain the observed corrugations? Our answer is yes, 
and there are two types of arguments that back this affir¬ 
mation: a) there is a non-negligible magnetic field in our 
Galaxy (Planck Collaboration et al. 2015 Beck 20151; and 


b) in the evolutionary models of the Galactic interstellar 
medium, the inclusion of the Galactic magnetic field is 
far from being irrelevant ([Martos fc Cox||1998 Santillan 


et al.|[l9M Gomez fc Cox|2004a|b|. 


5.1 Comparison between Observations and 
Models 

Our analysis establishes the existence of a vertical ve¬ 
locity field that shows a sinusoidal structure associated 
mainly with the observed spiral arms. How can the gen¬ 
eration of these structures be explained? Since the 1970s, 
different models have been proposed to explain the gene¬ 
sis of the spatial corrugations observed in our Galaxy and 


other external galaxies (see Alfaro & Efremov (19961 for a 
review of this topic to that date). But, regarding the vari¬ 
ations of the vertical velocity field, the first article that 


6 DIAGNOSTIC DIAGRAMS 

To complement the kinematic analysis, diagnostic dia¬ 
grams (DD) are included to study the ionization structure 
in the observed features. As it is well known, the ioniza¬ 
tion and excitation process can be determined by plots 
of [NH]/Ha versus [OHI]/H/?, depending on its location 
in such diagrams ( [Baldwin, Phillips fc Terlevich 19811. 
These ratios are the best options among all the possible 
emission line ratios (Veilleux & Osterbrock 19871, and 


therefore one of the most commonly and widely used DD. 

Before plotting the DD, the flux ratios are corrected 
for reddening following the relation 
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^Ha _ Fhol j^q0.4(Ajj^ — Ajjfl) _ Fho „ j^q0.4Av —fcg3) 

Ihp Fhp Fh(3 

(9) 

where I\ is the unreddened emission at a certain wave¬ 
length from the source, and F\ is the observed flux, af¬ 
fected by dust absorption. A\ = kxAv is the extinction 
factor at a particular waveband A, and kx and Ay are 
the extinction curve and the absorption in the V band, 
respectively. The Balmer decrement is a method for red¬ 
dening correction, deduced from the comparison of the 
observed Htr to H/? flux ratio with the theoretical value 


of 2.86 I Osterbrock 19891, based on case B recombination. 


Applying the Balmer decrement we have that. 


Av = 


log(2.86) - log{FHi3/FHa) 


0A{kHa — knp) 

So, the unreddened flux ratios can be obtained as 


( 10 ) 


log(^ 


1X2 1 


log(^ 

log(^ 


Fxi \ 
Fx2) 


-|- QAAvikxi — kx2) — 


( 11 ) 


where Ci = -fcgg ’ ~ log{2.86)Ci. 

The DD for each galaxy are found in Figures to 
together with the plots of AVz- Plotted individually, they 
have a better resolution than the above plots of AVz in 
Figures |3]-|^ To include more robust information, only 
those pixels with AVz less than the mean of the errors 
AVz are plotted. So they clearly show the relationship 
between the velocity and the Ha emission line flux peaks. 

The different slit position angles for each galaxy are 
represented with different symbols and grey colour scale, 
as indicated in the figures. The locations where the DD 
ratios are above the reference line for photoionization by 
young stars are shown in the panels with the AVz kine¬ 
matics. This could provide some clues as to whether the 
ionization mechanism is related to the presence of cor¬ 
rugations and the encounter of the density wave with a 
spiral arm. Within the galactic bore model framework 
( [Martos fc Cox|[l998| [Martos et al.|[T999P , we would ex¬ 
pect the presence of shocks at these locations. 

No concluding results are obtained from the DD, 
which show that photoionization by young massive stars 
seems to be the main ionization mechanism in most loca¬ 
tions. 

There is only a small portion of the gas that appears 
to be ionized by low-velocity shocks. In NGC 278 and in 
NGC 1058, at PA= 125°, this portion is more important. 
In NGC 278, for example, these low-velocity shocks would 
be located in the innermost central kpc, coinciding with 


Moreover, two different forms of behaviour appear in 
the [NII]/Ha vs. [OIII]/H/3 diagnostic diagrams. NGC 278 
and NGC 1058 have most of their pixels, for the three slit 
positions, concentrated at the same location in the DD, 
whereas for NGC 2500 and UGC 3574 different slit posi¬ 
tions occupy a different location in the DD, nearly cov¬ 
ering the theoretical curve separating the different ion¬ 
ization mechanisms. This behaviour may be related to a 
metallicity gradient, but a more detailed study, beyond 
the scope of this paper, would be needed. 

In summary, although the association between the 
gas ionization mechanisms and kinematic properties do 
not offer a conclusive pattern, some general characteris¬ 
tics are observed: a) most of the gas is ionized by high- 
energy photons; and b) in NGC 278 and NGC 1058 there 
is a not negligible number of points whose diagnostic dia¬ 
gram suggests ionization through low-velocity shocks. In 
NGC 278 these shocks seem particularly limited to the 
central region where some authors have observed signs 


of an inner merger (Knapen et al. 20041. Otherwise, we 


note the potential to explore this idea further with the 
integral field spectroscopic data sets that are becoming 


al. 


al. 


available, such as from the surveys CALIFA (Sanchez et 
2012| Garcfa-Benito et al.||2015 l, MaNGA (Bundy et 


2015p, and SAMI ( Croom|2012 1, and with much better 


resolution from data already becoming available from the 
MUSE instrument at the VLT (Bacon et al.|20i0|. 


the minor merger found by Knapen et al. (20041. 


7 SUMMARY 

In this work we extend the detection and study of kine¬ 
matic corrugations, already reported in the literature (e.g. 
Alfaro et al.|200T l, to a small sample of face on spirals. 

The galaxies studied here show a similar behaviour to 
NGC 5427 ( Alfaro et al.|[200T | at some locations, with a 
clear displacement between the velocity and the emission 
line peaks, where the approaching peaks of AVz occur 
in the convex side of the spiral arms, and the receding 
maxima are located behind the Ha flux maxima, in the 
concave side. This kinematical pattern is similar to the 
one expected in a galactic bore generated by the inter¬ 
action of a spiral density wave with a thick magnetized 
gaseous disk. 

In locations where AVz peaks are found to be related 
to Ha emission maxima, a strong velocity gradient cross¬ 
ing the arm is sometimes observed, as seen in Figures 
to[^ and described in Section 

This relation is not so clear in NGC 2500 and 
UGC 3574. However, in these galaxies the Ha emission 
is much fainter, and the Ha emission line could only be 
fitted in a limited number of points along the slit. This 
produces a bias, and no clear conclusions can be deduced 
for these two galaxies. 
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Figure 6. Top left: Diagnostic Diagram [NII]/Hq vs. [0III]/H/3 jVeilleux Osterbrock||1987[ l, for NGC 278. The resting plots 
are the detrended perpendicular velocity AV^ for the different slit positions. In this case, only those pixels with error 6Vz < ^Vz 
are plotted. The highlighted pixels correspond to those locations where the ionization mechanism is not photoionization. The Ha 
intensity is drawn in the right-hand ordinate, with a dashed line. A dotted line marking the zero value stands out as the global 
trend has been practically removed, with only the local oscillations of Vz remaining. The x-axis represents the distance to the 
reference pixel, in kpc units. 
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Figure 7. Analogous to Figure]^ but for NGC 2500. 


In NGC 278, at PA=143.5° there is a strong AVz 
peak in a weak Ha emission zone. This could be associated 
with corotation or with a recent merger in the inner disk. 
We have checked that the stronger velocity peaks do 
not necessarily coincide with the stronger emission line 
intensity peaks. 

From the analysis of the ionized emission line di¬ 


agnostic diagrams, we conclude that photoionization by 
young stars is the main mechanism, and only a small por¬ 
tion of the gas appear to be ionized by alternative mech¬ 
anisms, such as low velocity shocks. In this DD we see 
different behaviour for the two galaxies NGC 278 and 
NGC 1058 on one hand, and NGC 2500 and UGC 3574 
on the other hand. The former have most of the pixels 
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NGC 1058 



NGC 1058, PA = 125° 



galactocentric distance(kpc) 



galactocentric distance(kpc) 


Figure 8. Analogous to Figure]^ but for NGC 1058. 


concentrated at the same location of the DD, whereas for 
the latter they are distributed covering more of the pa¬ 
rameter space along the theoretical curve separating the 
different ionization mechanisms. This may be understood 
as a consequence of metallicity gradients, but further in¬ 
vestigations are warranted. 

The origin of kinematic corrugations is still a matter 


of debate. Corrugations are closely linked, as cause-effect, 
to the large scale star formation processes: density waves, 
tidal interactions, collisions of high velocity clouds with 
disk, or a galactic bore generated by the interaction of a 
spiral density wave with a thick gaseous disk, etc. Which 
mechanism is the origin of disk gas corrugations is still an 
open question. 
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Figure 9. Analogous to Figure]^ but for UGC 3574. 
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